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Abstract 



CN I The impact of the recent very precise CLAS and COMPASS gi/Fi data on polar- 

ized parton densities and higher twist effects is discussed. It is demonstrated that 
the low CLAS data improve essentially our knowledge of higher twist correc- 
' tions to the spin structure function gi, while the large COMPASS data influence 

mainly the strange quark and gluon polarizations. It is also shown that the un- 
certainties in the determination of the polarized parton densities are significantly 
i reduced. We find also that the present inclusive DIS data cannot rule out a negative 

polarized and changing in sign gluon densities. The present status of the proton 
spin sum rule is discussed. 



1 Introduction 



I One of the features of polarized DIS is that a lot of the present data are in the preasymp- 

O ■ totic region (Q^ ^ 1 - 5 GeV^ 4 GeV^ < < 10 GeV^). This is especially the case for 

the experiments performed at the Jefferson Laboratory. As was shown in [1], to confront 
^ , correctly the QCD predictions to the experimental data including the preasymptotic re- 

^ I gion, the non-perturbative higher twist (powers in 1/Q^) corrections to the nucleon spin 

structure functions have to be taken into account too. 

In this talk we discuss the impact of the recent very precise CLAS [2] and COMPASS 
[3] inclusive polarized DIS data on the determination of both the longitudinal polarized 
parton densities (PDFs) in the nucleon and the higher twist (HT) effects. These experi- 
ments give important information about the nucleon structure in quite different kinematic 
regions. While the CLAS data entirely belong to the preasymptotic region and as one 
can expect they should mainly influence the higher twist effects, the COMPASS data on 
the spin asymmetry Af are large data and they should affect mainly the polarized 
parton densities. In addition, due to COMPASS measurements we have for the first time 
accurate data at small x (0.004 < x < 0.015), which allow to determine the behavior of 
the PDFs at small x region and therefore to calculate more precisely the first moment of 
the nucleon spin structure gi. 
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2 NLO QCD analysis of the data 



The method used to extract simultaneously the polarized parton densities and higher 
twist corrections to the spin-dependent nucleon structure function gi is described in [1]. 
According to this method, the gi/ Fi and Qi/ Fi) data have been fitted using the 

experimental data for the unpolarized structure function Fi{x, Q"^) 

Fi{x,Q^) 

As usual, Fi is replaced by its expression in terms of the usually extracted from un- 
polarized DIS experiments F2 and R and the phenomenological parametrizations of the 
experimental data for F2{x, Q^) [4] and the ratio R{x, Q^) of the longitudinal to transverse 
cross-sections [5] are used. Note that such a procedure is equivalent to a fit to {gi)exp, 
but it is more precise than the fit to the gi data themselves actually presented by the 
experimental groups because here the gi data are extracted in the same way for all of the 
data sets. 

In Eq. (1) "LT" denotes the leading twist contribution to gi 

9i{x, Q\t = 9i{x, Q\qcd + h^'^ix, Q')/Q' + 0{M'/Q') , (2) 
where (5^)pqcd is the well known (logarithmic in Q"^) NLO pQCD contribution 

9.{x. Q^)pQCD = \ E elm + Ag) (1 + ^^SC,) + ^^AG ^], (3) 

and Q^) are the calculable kinematic target mass corrections [6], which effectively 

belong to the LT term. In Eq. (3), Aq{x , Q'^) , Aq{x , Q"^) and AG{x,Q'^) are quark, anti- 
quark and gluon polarized densities in the proton, which evolve in according to the 
spin-dependent NLO DGLAP equations. 5C{x)q^G are the NLO spin-dependent Wilson 
coefficient functions and the symbol ® denotes the usual convolution in Bjorken x space. 
Nf is the number of active ffavors (Nf = 3 in our analysis). h{x)/Q'^ in Eq. (1) corresponds 
to the first term in the (Aq^q/Q^)"^ expansion of higher twist contribution to gi. Its 
logarithmic dependence, which is not known in QCD, is neglected. Compared to the 
principal dependence it is expected to be small and the accuracy of the present data 
does not allow its determination. Therefore, the extracted from the data values of h{x) 
correspond to the mean for each x-bean. 

Let us discuss now how inclusion of the CLAS EGl proton and deuteron gi/Fi data 
[2] and the new COMPASS data on Af [3] influence our previous results [7] on polarized 
PDFs and higher twist obtained from the NLO QCD flt to the world data [8], before the 
CLAS and the latest COMPASS data were available. 

3 Impact of the new data on polarized PDFs and HT 

The new CLAS EGl/p,d data on gxjFx (633 experimental points) [2] and the recent 
COMPASS data on the longitudinal asymmetry A\ (15 experimental points) [3] are at 
very different kinematic regions. While the CLAS data are high-precision data at low 



^i(x,Q2)LT + /^,(,x)/g2 



(1) 



exp 
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Q^:{x - 0.1 - 0.6, ^ 1-5 GeV^, W > 2 GeV}, the COMPASS data are mainly 
at large (5^:{0.0046 < x < 0.57, ~ 1 — 55 GeV^} and the only precise data covering 
the low X region. Therefore, they will play a different role in the improvement of the 
determination of the polarized PDFs and higher twist effects. The new PDFs and HT and 
their uncertainties will be compared with those of LSS'05 determined from our previous 
analysis of the world data [8] available before the CLAS EGi/p, d and COMPASS'06 data 
have appeared. 

As the CLAS data are mainly low Q"^ data where 
the role of HT becomes important, they should help 
to fix better the higher twist effects. Indeed, due to 
the CLAS data, the determination of HT corrections \ 
to the proton and neutron spin structure functions, 
hP{x) and is significantly improved in the CLAS ^ 

X region, compared to the values of HT obtained from 
our LSS'05 analysis [7] in which a NLO(MS) QCD 
approximation for gi{x,Q^)ur was used. This effect 
is illustrated in Fig. 1. One can conclude now that 
the HT corrections for the proton target are definitely 
different from zero and negative in the x region: 0.1- 
0.4. Also, including the CLAS data in the analysis, 
the HT corrections for the neutron target are better 
determined in the X region: 0.2-0.4. Note that /i"(a;) at 
X ~ 0.5 was already fixed very precisely from the JLab 
Hall A data on the ratio Qi''^ / . We have found that 
the impact of the COMPASS'06 data on the values Figure 1. Effect of new data on the 
of higher twist corrections and their uncertainties is higher twist values, 
negligible. The only exception are the central values 
of HT at small x for both the proton and the neutron 

targets which arc slightly lower than the old ones. Note that this is the only region where 
the COMPASS DIS events arc at smaU Q^: 1-4 GeV^. 

The effect of the new data on the polarized PDFs and their uncertainties is demon- 
strated in Figures 2 and 3, respectively. The central values of both the {Au + Au) and 
{Ad + Ad) parton densities do not change in the experimental region (the correspond- 
ing LSS'06 curves can not be distinguished from those of LSS'05). As one can see from 
Fig. 2 the new data influence only the polarized gluon and strange quark sea densities 
(while the magnitude of strange sea decreases at x < 0.1, the gluon density increases at 
X > 0.1). As expected, the central values of the polarized PD are practically not affected 
by the CLAS data. This is a consequence of the fact that at low the deviation from 
logarithmic in pQCD behaviour of gi is accounted for by the higher twist term in gi 
in Eq. (1). So, the change of the central values of the polarized gluon and strange quark 
sea densities is entirely due to the new COMPASS data. On the contrary, the accuracy 
of the determination of polarized PDFs is essentially improved due to the CLAS data 
(the dashed curves in Fig. 3). This improvement is a consequence of the much better 
determination of higher twist contributions to the spin structure function gi, as discussed 
above. The impact of COMPASS data on the uncertainties for the PDFs is also shown in 
Fig. 3 (the solid curves). As seen, they help to improve in addition the accuracy of the 
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Figure 4. Expected uncertainties for NLO(MS) polarized PDFs after including the data set to 
be collected with CLAS12 experiment including statistical and systematic errors 

determination of the gluon and strange sea quark polarized densities at small x: x < 0.2 
for the gluons and a; < 0.1 for the strange sea. 

An essential further improvement (the dashed lines in Fig. 4) can be achieved after 
including in the analysis the data set to be collected with CLAS12 experiment [9] planned 
to be performed using a 12 GeV electron beam at Jefferson Laboratory, USA. 

At the end of this Section we would hke to mention that all results on the PDFs 
presented here have been obtained when 5 x-bins have been used to extract the HT 
values. Due to the good accuracy of the CLAS data, one can split the measured x region 
of the world data set into 7 bins instead of 5, as used up to now, and therefore, can 
determine more precisely the x-dependence of the HT corrections to gi. The numerical 
results of the best fit to the data using 7 x-bins are presented in [10]. It is important 
to emphasize that the central values for the PDFs(5 bins) and PDFs(7 bins) excepting 
the gluons are very close to each other. However, the uncertainties for the PDFs(5 bins) 
are smaller than those for PDFs(7 bins), especially for As(a;) and AG(a;). That is why 
we prefer to present here the PDFs and there uncertainties corresponding to 5 bins in x 
using for the HT values. 

4 The sign of the gluon polarization 

We have observed also that the present inclusive DIS data cannot rule out the solutions 
with negative and changing in sign gluon polarizations (see Fig 5a). The shape of the 
negative gluon density differs from that of positive one. In all the cases the magnitude of 
AG (the first moment of the gluon density) is small: [AGI < 0.4 and the corresponding 
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Figure 5. Strange quark sea densities xAs{x) corresponding to the fits with AG > 0, AG < 
and changing in sign a;AG 

polarized quark densities {Au + Au) and {Ad + Ad) are very close to each other. The 
corresponding strange sea densities are shown in Fig. 5b. Note, however, that the uncer- 
tainties for PDFs corresponding to the solution with AG < are larger than those in the 
case of AG > (for more details see [10]). In Fig. 6 the ratio AG{x)/G{x) calculated 
for the different AG(x) obtained in our analysis and using G(,x)mrst'02 [H] ioi the unpo- 
larized gluon density, is compared to the existing direct measurements of AG(a;)/G(a;). 
The error band correspond to statistic and systematic errors of AG(a;). The most precise 
value for AG/G, the COMPASS one, is well consistent with any of the polarized gluon 
densities determined in our analysis. One can see from Fig. 6 that in order to choose 
between gluons with positive and negative polarization direct measurements of AG{x) at 
large x : x > 0.3 are needed. 

5 The proton spin sum rule and spin puzzle 

Using the values for the singlet and gluon polarizations AE(Q^) and AG{Q'^) at Q'^ — 
1 GeV^ obtained in our analysis (MS scheme):AE = 0.207±0.039 and AG = 0.237±0.153 
we have found the following value for the spin of the proton ai Q'^ — 1 GeV"^: 

s.^l^ ^As(g2) + AG(g2) + L,(g2) = 0.34 ± 0.15 + l,{q'). (4) 

So, in order to satisfy the proton spin sum rule (4) the sum of the quark and gluon orbital 
angular momentum = + should be different from zero and positive. Note that 
the quark orbital momentum will be determined soon from the data using the forward 
extrapolation of the generahzed paron densities (GPD). 

Let us finally discuss the so called " spin puzzle" - the discrepancy between the values 
of the singlet polarization AS: 0.2-0.3 in the DIS region and 0.6 at low Q'^{Q^ ~ Aqcd) 
(see Fig. 7a). For better understanding of the situation it is useful to use the JET 
factorization scheme [12], in which AE(g^) does not depend on Q^. Then, in this scheme 
it is meaningful to directly interpret the singlet polarization AE as the contribution of 
the quark spins to the nucleon spin and to compare its values obtained in the DIS and 
low g^ regions. The value of AEjet obtained in our LSS'06 analysis of the DIS data is 
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0.26 ± 0.08. 

On the other hand the well known value of 0.6 
for AEiQ^ ~ A^cd) = + + ^^Isea is pr^ 
dieted in the relativistic constituent quark model 
(CQM) [13]. However, this model does NOT ac- 
count for the vacuum (quark sea) polarization. It 
was qualitatively shown in the instanton models 
[14, 15] that due to the non-perturbative vacuum 
spin effects the contribution of the sea quark po- 
larization to AS is negative. So, the value of AS 
in the non-perturbative region (Q^ 



^qcd) 



is 



really smaller than 0.6. Also, it was found from 
a combined analysis of forward scattering parity- 
violating elastic asymmetry data from G*^ and 
HAPPEx experiments at JLab, and elastic vp and 
vp scattering data from Experiment 734 at BNL, 
that the strange axial form factor G';|((5^), which 
is strongly related with As {G\{Q'^ = 0) = As), 
is negative in the region 0.4 < < 1 GeV"^ [16] 
(see Fig. 7b), i.e there is a strong indication that 
the strange quark contribution to AS at low is 
negative. In conclusion, we are very close to the 
solution of the so called "spin puzzle". 
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* COMPASS, high p, Q'>1GeV' 
T HERMES, (preiim.) 

■ SMC. high p,Q^>1GeV^ 




Figure 6. Comparison between the ex- 
perimental data and NLO(MS) curves 
for the gluon polarization AG{x)/G{x) 
at = 3 GeV^ corresponding to AG > 
0, AG < and an oscillating-in-sign 
xAG. 



Conclusion 

We have studied the impact of the CLAS and latest COMPASS data on the polarized 
parton densities and higher twist contributions. It was demonstrated that the inclusion 
of the low CLAS data in the NLO QCD analysis of the world DIS data improves 
essentially our knowledge of HT corrections to gi and does not affect the central values of 
PDFs, while the large COMPASS data influence mainly the strange quark and gluon 
polarizations, but practically do not change the HT corrections. The uncertainties in the 




Figure 7. A possible explanation of the nucleon's spin puzzle (a). Results of analysis for the 
strange axial form factor of the proton (b). 
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determination of polarized parton densities is significantly reduced due to both of the 
data sets. These results strongly support the QCD framework, in which the leading twist 
pQCD contribution is supplemented by higher twist terms of 0{K^qy>/Q'^)- 

Finally, one of the important messages coming from this analysis is that it is impossible 
to describe the very precise CLAS data if the HT corrections are not taken into account. 
Note that if the low data are not too accurate, it would be possible to describe them 
using only the leading twist term in gi (logarithmic in Q^), i.e. to mimic the power in Q"^ 
dependence of g\ with a logarithmic one (using different forms for the input PDFs and/or 
more free parameters associated with them) which was done in the analyses of another 
groups before the CLAS data were available. 
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